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Abstract 

We calculate the effect of the interaction between an optically active material and a Bose-Einstein 
condensate on the collective oscillations of the condensate. We provide explicit expressions for the 
frequency shift of the center of mass oscillation in terms of the potential generated by the substrate 
and of the density profile of the gas. The form of the potential is discussed in details and various 
regimes (van der Waals-London, Casimir-Polder and thermal regimes) are identified as a function 
of the distance of atoms from the surface. Numerical results for the frequency shifts are given 
for the case of a sapphire dielectric substrate interacting with a harmonically trapped condensate 
of ^^Rb atoms. We find that at distances of 4 — S/jm, where thermal effects become visible, the 
relative frequency shifts produced by the substrate are of the order 10~^ and hence accessible 
experimentally. The effects of non linearities due to the finite amplitude of the oscillation are 
explicitly discussed. Predictions are also given for the radial breathing mode. 
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I. INTRODUCTION 



The study of the force felt by an atom near a surface has recently become a popular subject 
of research (see for example and references therein) both in view of nanotechnological 

applications and because of the possibility of investigating fundamental forces at the 
submicron scale Q]. These studies extend previous investigations of the generalized van der 
Waals forces between macroscopic bodies ||5|] to more microscopic systems. Experiments 
with Bose-Einstein condensates near a surface 0, 0] are very promising in this respect 
because of the high accuracy of the measurements achievable with ultracold gases. 

Historically the interatomic forces were first introduced in 1873 by van der Waals in order 
to explain the deviations of the thermodynamic behaviour of real gases from the ideal laws. 
Only in 1930 London provided a first quantitative description of the force. He applied 
second order perturbation theory to the quantum Hamiltonian to evaluate the electrostatic 
interaction between two dipoles, induced by the fluctuation of the electromagnetic field. 
He succeeded in deriving the most famous power law l/(i^, usually known as the van der 
Waals-London interaction, where d is the interatomic distance. Almost twenty years later, 
in 1948, Casimir and Polder taking into account relativistic retardation effects within 
fourth order perturbation theory, found that, at large distances, the potential decays like 
rather than like 

The existence of these attractive forces between neutral atoms gives rise to analogous 
forces between a single atom and a dielectric or a metallic wall, as well as to forces between 
macroscopic bodies separated by small distances Q]. Between the end of the fifties and 
)eginning of sixties Lifshitz developed the general theory of the van der Waals forces 
13^ for inhomogeneous media, taking into account the quantum fluctuations of the 
electromagnetic field and their perturbation produced by the media. The theory describes 
the interaction between continuous media in terms of their dielectric functions and of the 
temperature. The van der Waals-London and Casimir-Polder forces are recovered as limiting 
cases of the more general theory, which also accounts for the effects induced by the thermal 
fluctuations of the electromagnetic field. Thermal effects provide the asymptotic behaviour 
of the interaction at distances much larger than the thermal wave length of photons. 

The purpose of the present work is to show that the study of the collective oscillations 
of a Bose-Einstein condensate of ultracold atoms can provide a useful probe of such forces. 
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Bose-Einstein condensates are very dilute samples and consequently the force produced by 
the surface can be safely calculated starting from the interaction felt by the individual 



atoms in the condensate 



ll|. In this paper we will be mainly concerned with the center of 



mass oscillation of a condensate confined by a harmonic trap. The center of mass motion 
has the main advantage of being independent of two-body interactions and its frequency 
consequently provides an ideal probe of the additional force generated by the wall. The 
effect of the force on the frequency of the center of mass oscillation depends on the distance 
of the condensate from the surface, on the temperature of the surface as well as on the 
optical properties of the surface and of the atoms in the condensate. The effect also depends 
on the thickness of the condensate which in actual configurations can be of the same order 
as the distance from the wall. In some limiting cases, corresponding to short and large 
distances, it is possible to obtain analytic or semi-analytic predictions that will be discussed 
systematically in our paper. They can serve as a natural guide for first quantitative estimates 
and as useful tests of more complete numerical approaches. 
The paper is organized as follows: 

In sect. 2 we derive some general results for the shift of the collective frequencies of a 
harmonically trapped Bose-Einstein condensate subject to an additional weak external force. 

In sect. 3 we employ Lifchitz theory to calculate the force acting on a dilute atomic gas 
placed near a surface. 

In sect. 4 we discuss the transition from the van der Waals-London potential, holding 
at very short distances from the surface and decaying like to the Casimir- Polder 

potential which accounts for relativistic retardation effects and holds at larger distances 
where it decays like 

In sect. 5 we discuss the transition from the l/d^ law characterizing the Casimir Polder 
potential to the behaviour of the potential at much larger distances from the surface (larger 
than the thermal wavelength of photons) where it decays like l/c?^ as a consequence of the 
thermal fluctuations of the electromagnetic field. 

In sect. 6 we discuss the role of the optical properties of the interacting media with speciflc 
emphasis to the case of the dielectric sapphire (^4/203) substrate and of the Rubidium atoms 
which form the condensate. 

In sect. 7 we provide results for the shifts of the collective frequencies as a function of 
the distance of the condensate from the surface, calculated at different temperatures. 
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Finally, in Appendix Al we provide results for the frequency shift of the radial compres- 
sion mode in very elongated harmonic traps. 



II. CENTER OF MASS OSCILLATION OF A TRAPPED GAS IN THE PRES- 
ENCE OF A WEAK PERTURBATION. 



After the first experimental realization of Bose-Einstein condensation in dilute vapours of 
alkali atoms the experimental and theoretical research in the field of ultacold quantum gases 
has grown in an impressive way (for general reviews see, for example, IJ, ll5|, |l6|, llll ) • Bose- 
Einstein condensed gases are dilute, ultracold samples characterized by unique properties 
of coherence and superfluidity. They give rise, among others, to interference phenomena 



m 



as well as to a variety of collective oscillations 
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22| . The relatively high 



density of Bose-Einstein condensed samples (compared to the one of non condensed trapped 
gases), as well as the possibility of imaging the atomic cloud after expansion, where the 
measured sizes become macroscopically large, permit to achieve very accurate determinations 
of the frequencies of the collective oscillations, thereby providing interesting opportunities 
for precision measurements. 

In this section we will calculate the effect of a weak perturbation on the center of mass 
oscillation of a harmonically trapped gas. The frequency of this oscillation, in the absence 
of the perturbation, coincides with the frequency of the trap. This result, also known as 
Kohn theorem 2^, is independent of the amplitude of the oscillation, interatomic forces, 
temperature as well as of the quantum nature of the sample. Deviations of the measured 
frequency from the oscillator value can then be used as a useful probe of additional inter- 
actions, like the one generated by the surface. In Fig.l we show a schematic representation 
of the geometry considered in the present work. A dilute gas trapped by the harmonic 
potential 



Vho{r) = -ujy + -uy + -uy (i) 

is placed close to an ideal wall generating a potential Vsurf{z) so that the full potential felt 
by each atom is given by 

V{v) = VU^) + Vsurf{z) . (2) 

In the following we will calculate the changes in the frequency of the center of mass oscillation 
along the z-th. direction caused by the presence of the surface. In the appendix we provide 



results also for the radial breathing mode in the case of very elongated axi-symmetric traps 

The equations for the center of mass coordinate Zcm = (l/^)(Ei=i^i) and the total 
momentum = (XliliP^i) along the z-th direction (A^ is the number of atoms of the 
trapped gas) can be written in the general form 

dZcM Pz 



dt Nm 

and 



(3) 



= -muolZcM - j c/m(r, (4) 



N dt 

i=l 

where we have used expression ((T]) for Vho and we have expressed the average of the external 
force dzVsurf in terms of the density n of the gas here assumed to be normalized to unity 
(J drn = 1). Notice that the equation for the total momentum is not affected by the two- 
body potential. In the absence of the perturbation {Vsurf = 0) one recovers the equations of 
the harmonic oscillator and the center of mass oscillation corresponds to a rigid oscillation 
of the density. The effect of the perturbation Vsurf can then be calculated, in first order 
perturbation theory, using the scaling ansatz n{r,t) = no{x,y,z — Zcuit)) in (jH), where hq 
is the equilibrium density of the gas. The equation for the center of mass coordinate finally 
takes the simple form 

d'^ZcM 2 " 1 



dt 



2 



-ulZcM ~~j drno{r)d^Vsurf{z + Zcm) (5) 



where, in the integral, we have made the substitution (z — Zcm) z. Equations and (jSj) 
have now a closed form. In the limit of small oscillations we can expand the perturbation 
term up to terms linear in Zcuit) and the equation for the center of mass takes again a 
harmonic form with the renormalized oscillator frequency 



ujI + ^ j no(r) dlVsurfiz) dr (6) 



It is worth stressing that result (jH)) is exact up to first order corrections in Vgurf and holds 
independently of the quantum or classical nature of the gas which instead determines the 
actual form of the density profile no. Eq.® can be further simplified by introducing the so 
called ID column density 

nl{z) = / dxdyno{x,y,z) (7) 



and using the fact that Vsurf depends only on the 2;-th coordinate. One then finds the result 

-hM = u.l + ^ P-oW%^^- (8) 
Equation (jHl) shows that the ingredients needed to calculate the effect of the perturbation 
on the shift of the center of mass frequency are the explicit form of the potential Vsurf, 
which will be discussed in the following sections, and the ID column density ((Tj). For a Bose- 
Eintsein condensate in its ground state the density profile is easily evaluated in the so called 
Thomas-Fermi approximation. This approximation holds for large condensates where, for 
gases interacting with positive scattering length, the kinetic energy of the trapped condensate 
can be neglected and one finds the analytic result [l^ no{x,y,z) = (/i — Vho{r)) when 
no > and elsewhere. Here /i = {huJho/ a / ahoY/^ is the chemical potential, is the 
number of atoms, uj^o = {tOx^y^zY^^ is the geometrical average of the three trap frequences, 
ttho = ^/^/ifTT-^ho) is the oscillator length and g = Anh'^a/m is the interaction coupling 
constant fixed by the s-wave scattering length a. Integration of n with respect to x and y 
yields 

15 1 / ^2X2 



where is the Thomas-Fermi radius in the z-th direction, fixed by the relation /i = 
mujlRl/2. Typical values for Rz in standard experimental conditions can be of the order of 
a few microns j^]. 

Equations (jH)) and are the key result of this section and will be used in the second part 
of the paper to calculate the effect of the force produced by the surface on the frequency shift 
of the center of mass oscillation. The integral (jH)) should be in general carried out numerically. 
There are, however, important cases where the integration can be found analytically. This 
is the case, for example, if one can approximate the potential with the algebraic decay 
(z) = —an/ (d + z)" where d is the distance between the wall and the center of the 
harmonic trap j2J]. In this case we find the result 

2 2 20 04 1 , . 

'^^^ = ''^-^ (1^ 

and 

2 _ 2 15 as [10 - 6 X - 3(1 - x'? In ((1 - + x))] 

^CM-^z 4 ^^5 {I-X^YX" ■ ^ ^ 

for the most relevant n = 4 and n = 3 cases, respectively, where x = Rz/d. Results 
(fTU|) and (fTT|) explicitly point out the role played by the finite size of the condensate. In 



particular, only in the limit i?^ ^ (i (x <^ 1) will the above equations approach the power 
law behaviours and for the frequency shift expected for a point-like condensate. 

The above results have been derived in the linear regime. If the amplitude of the oscilla- 
tion is comparable to the distance from the surface non linear effects become important and 
modify the value of the frequency shift. Non linear effects result in a frequency shift as well 
as in the occurrence of higher harmonics. They can be calculated starting from the general 
equations Q and for the center of mass coordinate. Notice that, due to the harmonic 
nature of V/io, both effects are absent if Vgurf = 0. Due to the perturbative nature of the 
problem we can use the periodic law Zch[ = acosiucMi) to evaluate the integral of ©, 
where a is the amplitude of the oscillation. One can then introduce the periodic function 

Q{t) = — I drno{r)d^Vsurf{z + acos{uJcMt))dz (12) 
ni J 

which can be conveniently expanded in a Fourier series. For the calculation of the frequency 
shift only the term proportional to cos{uJcMt) is important and one can write 

2it/lucm 

Q{t) cos{uJcMt)dt. (13) 

Expansion of pHj) up to terms quadratic in a and use of the ID column density finally yield 
the result 

- = - nl(z)^^dz + ^ J_^^ nl(z)^^dz (14) 
for the frequency shift which generalizes eq.(jHl) by including the non linear correction in a?. 

III. FORCE BETWEEN THE SURFACE AND A SINGLE ATOM. 

The interaction produced by a surface on a single atom includes, in its general form, non 
trivial relativistic effects as well as quantum and thermal fluctuations of the field. The force 
can be calculated starting from the general theory devoloped by Lifchitz for the free energy 
associated with inhomogeneous media. For two plates separated by a distance d and in 
thermal equilibrium with the thermal radiation at temperature T, the force Fgurf per unit 



2 2 _ ^CM I 

vra /n 
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surface can be written as 



surf 



(d) 



16 TT d? 



V 



n=l 



(£10 + 1) (£20 + 1) ^ 

(£10 - l)(e2o - 1)' 
exp 



-1 



dx+ 



+ 



(15) 



+ 



{si + p ei){s2 + p 62) f2pind 

exp 



dp 



[si - p ei){s2- p 62) ' V c 
where, with respect to the usual presentation of the force, we have separated a first term 
behaving hke ksT/d^ from the other terms with n > 1. In eq. ()15|l we have defined 



Sl 



- 1 + p2 ^ S2 = ^£2 - 1 + 



(16) 



where £1 = ei{i ^„) and £2 = ^2(^ ■Cn) are the relevant dielectric functions of the two plates, 
evaluated for imaginary values u = i C,n of the frequency, with ^„ = 2iTkBTn/h. The 
function e{i ^) is a monotonic and real function of ^, obtained by analytic continuation of 
the non monotonic and complex function e{uj) evaluated for real values of the frequency 
u. Furthermore in eq.()15|) we have introduced the static values of the dielectric functions 
£^10 = ^i(O) and £20 = £^2(0). Some general properties of these functions will be illustrated in 
sect VI. 

The force produced by a surface on a single atom can be easily derived by assuming that 
one of the two walls (wall 2) is made of a very dilute material (for example a gas) so that 
one can expand eq. ljlSp for small values of 62 — 1 = Anna where n is the density of the gas 



and a is the corresponding atomic polarizability. The potential Knr/ felt by a single atom 

1 the forr 

2p^nd'' 



is then simply given by Vsurf = 



-Fsurf I n and can be written in the form 



a a 



ysurf{d) 



KT 



Ad^ (£0 + 1) - — , ji 
where we have introduced the function f{p) defined by 

(^1 



exp 



f{p) dp. (17) 



(18) 



[si +p) - {si+pe] 

In eqs. (jl7|) and (|18p we have omitted the index 1 in the dielectric function characterizing 
the optical properties of the surface while and are the static values of the dielectric 
function and of the atomic polarizability, respectively. Equation (fTTj) permits to describe 
both the interaction generated by a dielectric and by an ideal metal (in the latter case the 
function e(a;) behaves like l/uj at small ui and (eo — l)/(eo + 1) = 1- 
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Eqs. fjl5|) and p7|) have been derived assuming that the whole system is in thermody- 
namic equihbrium and in particular that its components have the same temperature T. In 
the experiments with Bose-Einstein condensates the atomic gas is cooled down to extremely 
low temperature, while the substrate remains at room temperature. Result (jl7j) is still ex- 
pected to be valid in this case, T being the temperature of the substrate. In fact atoms, 
independently of being at zero or at room temperature, do not contribute to the thermal 
radiation, their lowest excitation energies being much higher than ksT. 

It is useful to identify the relevant length scales of the problem (see fig. 2). A first impor- 
tant length is the thermal wavelength of the photon: 

A. ^ ^ (19) 

where ks is the Boltzmann constant. At room temperature (T = 300i^') one has At = 7.6fim. 
For distances d larger than At the force is dominated by the thermal fluctuations of the 
electromagnetic field. In this limit the leading contribution is given by the first term in 
eq.()17|). due to the occurrence of the exponential factor in the other terms, and the potential 
takes the characteristic form 

T . . ks T ao {so - 1) 

Notice that in this regime only the static value of the dielectric and polarizability functions 
contribute to the force. It is also worth noticing that in this limit the force is independent 
of the Planck constant as well as of the velocity of light. 

For distances d smaller (or of the same order) than At the quantum fluctuations of the 
electromagnetic field become important and the terms of the sum (fTTj) should be taken into 
account. If (i <^ At and T is sufficiently small, one can replace the sum with an integral 
(see sect. IV). In this regime one can introduce a second length scale, hereafter called Xopt- 
This length is related to the relevant wavelengths characterizing the optical properties of 
the interacting media. One can identify two distinct regimes. For distances d <^ Xopt 
the potential exhibits the familiar l/# van der Waals-London dependence. In the interval 
Xopt d <^ Xt one instead enters the Casimir-Polder regime where the potential decays 
like 1/d^. The possibility of identifying the Casimir-Polder regime depends crucially on the 
value of the temperature. The temperature should be in fact sufficiently low in order to 
guarantee the condition At ^ Xopt (see discussion at the end of sect. V). Differently from 



the thermal wawelength, given by the simple expression At, the explicit evaluation of Xopt 
is not immediate because the optical properties of the media are not characterized by a 
single frequency (see sect. VI ). In the case of Rubidium atoms interacting with a dielectric 
sapphire substrate we find Xopt ~ O.l^m. 



IV. THE FORCE AT SHORT DISTANCES: FROM THE VAN DER WAALS- 
LONDON TO THE CASIMIR-POLDER REGIME. 

As anticipated in the preceding section at distances smaller than the thermal wavelength 
also the terms in the sum ()17|) contribute to the force and, if <^ Ay and the temperature 
is sufficiently low, the sum is conveniently replaced by an integral: J2n ^ {h/2nkBT) f d^. 
One then finds the useful expression 

- -s^^if^^r* " - ^) pi) 

for the potential, also called short range (SR) approximation. It is worth noticing that result 
(12 1|) does not depend on the value of T and coincides with the zero temperature limit 
of the general form (fT7j) of the potential. Equation (PT|) admits two important limits. For 
large d one finds the most famous expression 



also know as the Casimir-Polder law, where we have introduced the function l28l 



(e + 1) f 1 4 - (g + l)^e Arsh^/(^^ [l + e + 2e{e- if] 

(£-l)l3 2(5-1) 2(5-1)3/2 

+ — 1 ^ (23) 



which is equal to 1 in the case of an ideal metal (^o = oo). For sapphire (eo = 9.4) one finds 
0(£o) = 0.8 

In the opposit limit of small d one instead finds the different law 



also know as the van der Waals-London (VL) interaction. The coefficient characterizing this 
force is fixed by the optical properties of the media. The matching between the two laws 
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and provides a useful estimate of the optical length Xopt according to 
3 c ao (f^o ~ 1 



Xopt = — 0(go) 

/ £0 + J- 



£i!iL^ d^l \ (25) 



£(zO + 1 

In fig. ISlwe show the potential (pT|l (dashed-dot line) calulated using the dielectric functions 
of Rubidium atom and of sapphire. The asymptotic laws V^^j- = —Q8.1/d^[nK{fim)^] for 
the van der Waals (see eq. p^ ) and V^^j = — 8.34/ d'^[nK{^m)'^] for the Casimir- Polder (see 
eq. (1221)) potentials, as well as the asymptotic thermal law V^^j = — 2.86/ d^[nK{fim)^] 
evaluated at T = 300K (see eq.fOUj)). are also shown {d is here expressed in microns). It is 
worth noticing that the potential (j?T|) approaches the asymptotic Casimir-Polder force only 
at distances significantly larger than Xopt- 

V. THE FORCE AT LARGE DISTANCES AND THE STATIC APPROXIMA- 
TION: QUANTUM VS THERMAL EFFECTS. 

At large distances from the wall the potential ()17|) can be safely evaluated in the so 
called static approximation (SA) which consists of replacing the dielectric functions e and 
a with their static values Eq and ao respectively. In the static approximation the potential 
generated by the surface can be written in the useful form 



V^fi") - -'^ Oiv) (26) 



where we have introduced the relevant scaling variable 



y = f (27) 

AT 



and the function 



G{y) 



.3 (^0 + 1) r ..3 



(28) 



In eq.((2HI) /o(p) is the static limit of eq. lfTRjl . obtained by replacing e with Eq, while the 
function g is defined by 

e"" [1 + 4 e"" + 

lb a = 

n=l 

The asymptotic behaviours of the function G{y) are well established: for ?/ — > one finds 
G{y) —>■ 3(j){eo)/{2 n y). Viceversa for y ^ 00 one finds G{y) — > 1. For a fixed value of T the 
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two limits apply, respectively, to distances larger and smaller than the thermal wavelength 
(fTUj) . It then follows that the static approximation (j^Bj) provides the transition from the 
Casimir-Polder law (j221) holding for Xopt ^ d <^ A-p, to the thermal law ^2U^ holding for 
d ^ At- In fig. 4 we show the scaling function G{y) together with the two asymptotic limits 
?/ — s> and y — s> oo in the case of sapphire where we have used the value cq = 9.4. The 
function G{y) depends rather smoothly on the actual value of the static dielectric function, 
provided the value Eq is significantly larger than 1. 

In fig. 3 we compare the prediction for the potential given by the static approximation 
with the exact result evaluated starting from eq. p7|) . The comparison shows that for dis- 
tances larger than ~ 2/im the static approximation is rather accurate. This explains why 
the predictions of this approximation for the frequency shifts are also accurate (see figs El 
andEI). 

At shorter distances the static approximation fails and one should instead use the short 
range description developed in sect. IV. If At ^ Xopt the long range description based on 
the static approximation matches the short range description of sect IV in the Casimir- 
Polder regime ()22|1 . In this case the matching procedure completes the determination of the 
potential at all distances. If instead the condition Ay ^ Xopt is not well satisfied one should 
calculate explicitly the potential using the general expression ()17|) and the potential will not 
exhibit the intermediate Casimir-Polder behaviour ()22|) . For example in the case of fig. 3, 
where the potential is calculated at T = 3007^, the matching between the two curves is not 
very good and consequently the system never exhibits the Casimir-Polder behaviour j^^]. 

VI. OPTICAL PROPERTIES OF THE ATOM AND OF THE SUBSTRATE 

The functions e{i^) and a{i^), which are crucial ingredients for the calculation of the 
potential generated by the surface, obey the Kramers- Kronig relations 



where e" is the imaginary part of the dielectric function e = e' + ie" evaluated on the real 
axis. For the polarizability one has the analogous relationship 




(30) 




(31) 
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where a" is the imaginary part of the atomic polarizabihty a = a' + ia" . 

The imaginary part of the dielectric function can be related to the measurable real 
an imaginary parts of the refraction index n{ijj) = n'{uj) + in"{uj) according the equation 
e"{uj) = 2n' {uj)n" {uj) . The values of €"{u!) (data taken from [s^) and e(^0 sapphyire are 
reported on figs. [7| and |S1 From the optical point of view, crystalline sapphire is an uni- 
axial crystal and one should distinguish between the ordinary and extraordinary refraction 
indices. Since the optical properties for sapphire are not completely known, in our calcu- 
lations we have used only the data for ordinary waves. We expect that this approximation 
will not significantly affect the analysis in both crystalline and melted sapphires. In fig. IHl 
we report the corresponding values for a(i^) for ®^Rb js^. The static value of the polar- 
izabihty is Oq = 47.3 X 10~^^cm^. Figure [7| shows that the optical properties of sapphire 
are rather complex so that the actual value of Xopt governing the transition from the van 
der Waals-London to the Casimir-Polder regimes cannot be simply inferred from the form 
of e", but requires the explicit calculation of the integral (j25|) . The typical "two plateau" 
behaviour exhibited by e{i^) in solid dielectrics is the consequence of the concentration of 
the strength e"{uj) in two distinct regions of frequencies (see figs [3 and |H1) . This behaviour 
is absent in the atomic polarizabihty (see fig. Using eq. it is possible to calculate 
the value of Xopt- In the case of the sapphire substrate interacting with Rubidium atoms we 
find Xopt = O.lyum. 



VII. RESULTS FOR THE SHIFT OF THE CENTER OF MASS OSCILLATION 

We are now ready to calculate the shifts of the center of mass frequencies of the condensate 
caused by the force generated by the surface. We will consider a condensate of finite Thomas- 
Fermi radius located ad distance d from the surface (see fig.l). We will evaluate the 
fractional frequency shift 

7 = ^^1^^ , (32) 

with uJcM given by eq.®, as a function of the distance d at different temperatures . As 
discussed in the preceding sections, the behaviour of the potential Vsurf at short distances 
is temperature independent and coincides with the zero temperature value. Conversely, 
at distances of the order of the thermal wavelength, or larger, the potential exhibits an 
important temperature dependence. This behaviour is reflected in the shift of the center 
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of mass frequency of the condensate for which we provide our predictions in fig. 10 at 
different temperatures. The curves of this figure have been calculating using the general 
espression (II 7|) for the potential. In the calculation we have used the values c<j/27r = 220Hz 
for the harmonic trapping frequency along the z-th direction and = 2.5fim for the 2;th 
radius of the condensate. These are typical values employed in current experiments with 
Bose-Einstein condensates. We have used a sapphire substrate and the condensate is made 
of ^^Rb atoms. Our predictions show that at distances ~ 4 — 8fim, where the thermal 
effects start becoming important at room temperature, the frequency shifts are of the order 
of 10~"^. These results are promising in view of the possibility of systematic measurements of 
the interplay between quantum and thermal effects in the generalized van der Waals forces. 

It is also interesting to compare the predictions for the shifts of the collective frequencies 
with the results obtainined using the static approximation, holding at large distances (see 
sect.V). The comparison (see figs. ElandlHl) shows that the static approximation provides 
an accurate decription of the shifts for distances larger than 4/im, even at zero temperature. 
At shorter distances the deviations become more important. One should actually notice 
that at short distances the width of the condensate plays an important role amplifying the 
effects of the short range component of the potential. 

Let us finally discuss the effects of non linearities on the frequency shift. If the amplitude 
of the oscillation is not small additional corrections due to the external force should be taken 
into account (see sec. II). The corresponding corrections can be estimated using result (|TH) . 
For example, by choosing a = 0.5/im and making the same choices for the other parameters 
{Rz = 2.5 and uJz/2tt = 220Hz), we find that the prediction for the frequency shift ()H2|1 is 
increased, in absolute value, by ~ 20% at d = 6/im. This result shows that for larger choices 
of the amplitude of the oscillation and for smaller values of d the effects of non linearity can 
become very important. In general they should be calculated starting directly from eq. ()12|) 
and (Uni). 

VIII. CONCLUSIONS 

In this work we have applied Lifshitz theory of generalized van der Waals interactions to 
investigate the effect of a substrate on the collective oscillations of a trapped quantum gas. 
We have first developed the general theory for the frequency shifts of the center of mass 
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oscillation caused by a perturbative force of general form. We have hence characterized 
the potential generated by the surface identifying various relevant regimes. A first approx- 
imation, holding at short distances d from the surface, permits to describe the transition 
from the van der Waals-London law to the Casimir- Polder l/d"^ regime characterized 
by quantum and relativistic effects. A second approximation, holding at larger distances, 
permits to describe the transition from the Casimir- Polder law to the large distance ksT/d^ 
behaviour, dominated by the thermal fluctuations of the electromagnetic field. We have also 
identified the relevant length scales of the problem and explored the matching conditions 
for the various approximations. Due to the finite size of the condensate the calculation of 
the frequency shifts requires a numerical integration, although analytical expressions can 
be obtained in some limiting cases. The calculation requires the explicit knowledge of the 
optical properties of the interacting media and in particular the dielctric function of the 
susbstrate and the atomic polarizability of the atoms in the condensate. We have also ex- 
ploited the effects of nonlinearity and shown that, due to the finite width of the condensate, 
they can provide large corrections to the frequency shifts and should be consequently taken 
into account in the comparison with future experiments. We have finally derived (see the 
Appendix) results for the frequency shifts of the radial breathing mode. 

Our predictions are rather promising and suggest that Bose-Einstein condensates can 
become efficient sensors of very weak forces. In particular the oscillation of the condensate 
permit to explore the interplay between quantum and thermal fluctuations of the electro- 
magnetic field in the relevant region of 4 — 8//m. For such distances the relative shifts are 
predicted to be of the order of 10~^ and are hence measurable experimentally. 
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X. APPENDIX: FREQUENCY SHIFT OF THE RADIAL BREATHING MODE. 



The radial breathing oscillation in elongated traps is a fundamental mode exhibited 
by Bose-Einstein condensates. For an interacting Bose-Einstein condensate, harmonically 
trapped in two dimensions, this mode exactly occurs at twice the trapping radial frequency 
(we assume here symmetric trapping: uo^ = uo^ = uj±)., irrespective of the amplitude of the 
oscillation, number of atoms and value of the scattering length j^. In a 3D configuration 
this mode can be realized in the limit of a very elongated axi-symmetric trap {uoy <^ uj±) 
corresponding, in the geometry of fig. 1, to a condensate very elongated along the direc- 
tion orthogonal to the figure. The radial breathing oscillation has been measured with high 



accuracy m 
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The radial breathing mode can be described by deriving the time-dependent equations 
for the average square radius (x^ + z^) within Gross-Pitaevskii theory 17[ in 2D. After some 
straightforward algebra one finds the equation 

+ z') = --(r • drV) + —,{pI+pI) + -{V,^,o,y) (33) 

where V2-body is the average values of the mean field interaction energy and V is the external 
potential. By introducing the energy per particle E/N = {V) + {p^ + pi) /2m + {V2-body), 
which is conserved in time, and separating in V the axi-symmetric harmonic potential Vho 
from the surface potential Kur/, eq.(jSSl) can be rewritten as: 

^(^2 + z^) = ±E- Aul{x' + Z') - -{V^urf) - -^zd^Vsurf) (34) 

ctt^ m m m 

In the absence of the surface perturbation (Kwr/ = 0) eq. ljH^ gives oscillating solutions with 
ijj = 2u±. 

The effect of the perturbation is simply calculated noticing that the unperturbed so- 
lution corresponds to a scaling transformation where the density varies as n{x, y, t) = 
a^riQ^x / a,y / a) and no is the equilibrium density distribution. The time dependent scal- 
ing parameter a fixes the value of the square radius according to (x^ + z^) = a'^{x^ + z^)o- 
Using the scaling transformation to evaluate the integrals (Kur/) = / drn{x,y,t)Vsurf{z) 
and {zdzVsurf) = / dvn{x,y^t)zdzVsurf{.z) one can rewrite eq. ljMj) as a simple equation for 
the parameter a. In the following we will be interested in the limit of small amplitude 
oscillations. By expanding the integrals up to first order in (a — 1), one finally obtains the 
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result 



dznlzdzzdzVsurf + 2 / dznlzdzVsurf (35) 
■Rz J~Rz 



qZUz V surf 



m{x'^ + Z^) 

for the frequency shift of the radial breathing (B) mode, where tzq is the ID column density 
([7j). Similarly to the case of the shift of the center of mass frequency (see eq. flS^ ). also result 
fl35|) is exact up to first order corrections in the perturbation. In fact the radial scaling ansatz 
is an exact solution of the Gross-Pitaevskii equation in 2D. It is also a good approximation 
in 3D if the harmonic trap is sufficiently elongated. Result ()35|1 for the radial breathing 
mode may provide a further tool to investigate the effects of the generalized van der Waals 
force on a Bose-Einstein condensate. 
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Substrate 





FIG. 1: Schematic figure of the substrate-trapped gas system. Gravity is oriented along z. 
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FIG. 2: Relevant length scales Xopt and At caracterizing the generalized van der Waals potential Vsurf- 
The power law behaviour of Vsurf (d) is shown in the appropriate ranges of validity as a function of the 

distance from the surface. 
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FIG. 3: The atom-surface potential is shown using the exact formula of ea. (|17|l (solid line), the short range 
approximation H21|) (dash-dotted line) and the static approximation (|26|l (dashed line) . The asymptotic van 
der Waals-London (~ l/d^), Casimir-Polder (^ 1/c?'*) and high T (~ l/rf"*) potentials are also shown (dotted 
lines). The curves were obtained for a sapphire substrate at 300/^ and for '^^Rb atoms in the condensate. 



21 
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FIG. 4: The function G{y) (solid line) is shown as a function of the scaling variable y = d/Xx together with 
its asymptotic behaviours (dotted lines). The function was calculated for the sapphire substrate (see text). 
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FIG. 5: Fractional frequency shift (|32|l of the center of mass oscillation (solid line) as a function of the 
distance from the surface, compared with the prediction of the static approximation (dashed line). The curves 
were obtained for a sapphire substrate at 300K and for ^^Rb atoms in the condensate. The unperturbed 
frequency was uJz/2t: ~ 220Hz while the radius of the condensate was Rz — 2.5fim. 
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FIG. 6: Same as in fig.5, at T = OK. 
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FIG. 7: Immaginary part of the sapphire dielectric function e"{oj) on the real axis of frequences. 
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FIG. 8: Sapphire dielectric function s{i ^) evaluated on the imaginary axis. 
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FIG. 10: Relative frequency shift H32|) of the center of mass oscillation calculated at T = OK (dash-dotted 
line), T = 300 A' (solid line), T = 400 AT (dotted hne) and at T = 600 A' (dashed hne). The curves were 
obtained for a sapphire substrate and for ^^Rb atoms in the condensate. The unperturbed frequency was 
uJz/2tt = 220Hz while the radius of the condensate was Rz = 2.5/im. 
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